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This  report  presents  a  description  of  a  conputer  program  to  be  operated 
interactively  in  the  INTERCCM  environment  on  the  CDC  GfBER  coirputer  system. 
This  program  is  the  interactive  graphics  version  of  the  CRKGRO  program  devel¬ 
oped  under  a  research  program  entitle  "Improved  Methods  for  predicting 
Spectrum  Loading  Effects."  This  program  was  administered  by  the  Flight 
Dy-namics  Laboratory  of  the  Air  Force  Wright  Aeronautical  Laboratory,  Wright- 
Patterson  Air  Force  Base,  Ohio,  under  Contract  F33615-77-C-3121,  Project 
2401,  "Structural  Mechanics,"  Task  240101,  "Structural  Integrity  for  Military 
Aerospace  Vehicles,"  Work  Unit  2401020,  "Inproved  Methods  for  Predicting 
Spectrum  Loading  Effects." 

This  research  program  was  primarily  conducted  by  personnel  from  the 
Fatigue  and  Fracture  Mechanics  Group,  Dynamics  Technology,  Structure  Systems, 
supervised  by  George  E.  Fitch,  Jr.,  supervisor,  Joseph  S.  Rosenthal,  manager, 
and  Dr.  Leslie  M.  Lackman,  director.  James  B.  Chang  was  the  program  manager 
and  principal  investigator.  R.  Fink  was  responsible  for  the  initial  program 
development . 
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Section  I 


ir^TRODUCTION 


The  INTERACTIVE -CRKGRO  program  is  a  streamlined  fatigue  crack-growth 
analysis  conputer  code  that  operates  interactively  in  the  INTERCC^lfl^ 
environment  on  the  CDC  CYBER  conputer  system.  This  orogram  is  the  inter- 
active  version  of  the  detailed- level  crack- growth  analysis  program,  CRKGRO^ 
wtih  an  additional  capability  of  iiputting  user-coded  stress  intensity  factor 
solutions  used  for  the  crack-growth  analysis. 

In  recent  years,  crack-growth  analysis  has  become  inportant  throughout 
the  life  cycle  of  any  aircraft  system  since  the  Air  Force  philosophy  for  achie¬ 
ving  the  structural  integrity  of  airplanes  has  enployed  damage -tolerance  and 
durability  requirements,  A  crack-growth  program  with  an  interactive  capabi¬ 
lity  will  provide  realtime  feedback  schemes  for  users  to  observe  and  influence 
change  on  the  particular  growth  history,  and  will  certainly  enhance  evaluation 
of  crack-growth  behavior  in  any  phases  of  design,  even  as  early  as  the  pre¬ 
liminary  design  stage. 

The  crack-growth  methodology  implemented  in  CRKGRO  is  based  on  linear 
elastic  fracture  mechanics  principles  with  the  assertion  that  the  primary 
parameter  which  controls  the  growth  of  a  crack  is  the  stress  intensity 
factor  range  (AK)  at  the  crack  tip.  Although  there  is  a  wide  variety  of  stress 
intensity  factor  (K)  solutions  available  in  the  program  for  various  crack 
configurations,  the  option  to  enable  the  user  to  input  the  user-coded  stress 
intensity  factor  solution  will  undoubtedly  add  to  the  program  a  great 
versatility,  especially  in  the  evaluation  of  a  new  stress  intensity  factor 
solution  during  the  K-solution  development  stage. 

This  manual  presents  information  needed  to  use  INTERACTIVE -CRKGRO  and 
provides  input  and  output  guidance  to  the  users.  Users  are  advised  to  refer 
to  CRKGRO  User's  Manual f 2)  for  more  detailed  background  of  crack  growth 
analysis  procedures  and  input/output  descriptions. 
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Section  II 


C(M>UTER  PROGRAM  DESCRIPTION 


INTERACT I VE-CRKGRO  uses  the  CDC  time-sharing  software  package  INTERCOM 
and  the  TEKRONIX  plotting  software  package  PLOT- 10  f^)  to  provide  an  inter¬ 
active  capability  of  the  crack-growth  program,  CRKGRO.  In  addition,  the  CDC 
program  library  management  software  package  UPDATE  (4)  provides  users  access 
to  the  INTERACTIVE-CRKGRO  source  program  library  to  implement  an  additional 
feature  of  allowing  users  to  insert  user-defined  stress  intensity  factor 
solutions  into  the  INTERACTIVE-CRKGRO  subroutine  for  subsequent  conputation. 
INTERACTI\'E-CRKGRD,  contains  the  following  features  not  available  in 
CRKGRO: 

1.  Allows  inputs  from  a  specified  data  set  and  teminal  instead  of 
from  punched  cards 

2.  Allows  tabular  outputs  to  be  displayed  on  the  terminal  screen  and 
provides  hard  copies 

3.  Allows  graphical  outputs  to  be  displayed  on  the  terminal  screen 
and  provides  hard  copies 

4.  Allows  users  to  input  user-defined  stress  intensity  factor 
solutions  in  a  self-contained  format  on  site 

INTERACTIVE-CRKGRO  is  written  in  FORTRAN  IV  and  programmed  in  modular 
form  with  the  placement  and  retrieval  of  data  through  common  statements 
and  local  files.  The  main  modules  consist  of  nine  subroutines:  INPUT 
OUTPUT,  CRIT,  aCCNT,  GROW,  PTPAPM,  PLOT,  KIOXX,  and  K20XX.  A  brief 
description  of  each  of  these  modules  follows: 

INPirr  Reads  the  input  data  consisting  of  crack  growth  rate  constants, 

and  load  interaction  parameters,  fracture  properties  and  load 

or  stress  spectrum,  etc 

OUTPUT  Prints  the  echo  of  the  input  data 

CRIT  Computes  the  critical  crack  length 
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GfCCNT  Range-pair  counts  the  stress  spectrum 

GROW  Conputes  the  crack-growth  history  data 

PTPARM  Plots  the  crack-growth  history  data  in  four  forms 

PLOT  Plots  the  grid,  scaling,  and  data  point  for  PTPARM 

KlOxx  Stress  intensity  factor  library  for  part-through  cracks 

K20xx  Stress  intensity  factor  library  for  through-the-thickness  cracks 

A  list  of  the  code  number  assigned  to  each  crack  configuration  contained 
in  the  present  K- library  is  shown  in  Figure  1. 
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CODE  NO. 


DESCRIPTION 


SURFACE  CRACK,  CENTERED 


ONE  CORNER  CRACK  FROM 
CENTERED  OPEN  HOLE 


TWO  CORNER  CRACKS  FROM 
CENTERED  OPEN  HOLE 


ONE  CORNER  EDGE  CRACK 


USER-DEFINED  PART-THROUGH- 
CRACK  STRESS  INTENSITY  FACTOR 
SOLUTION 


THROUGH-CRACK,  CENTERED 


ONE  THROUGH-CRACK  FROM 
CENTERED  OPEN  HOLE 


TWO  THROUGH-CRACKS  FROM 
CENTERED  OPEN  HOLE 


ONE  THROUGH-EDGE  CRACK 

TWO  THROUGH-EDGE  CRACKS 

ASTM  STANDARD  CT  SPECIMEN 

USER-DEFINED  THROUGH-CRACK 
STRESS  INTENSITY  FACTOR  SOLUTION 


GEOMETRY 


Figure  1.  Crack  Library 


Section  III 


INPUT  AND  (XITPUT  DESCRIPTION 


INTERACT I VE-CRKGRO  accepts  input  data  from  two  sources.  The  primary 
source  is  the  terminal  keyboard,  through  idiich  the  user  types  in  data  in 
response  to  each  input  query  from  the  program.  The  alternate  form  of  input 
is  through  an  external  data  file.  This  file  contains  the  card  images  of  the 
input  data  that  would  normally  be  supplied  via  the  terminal  keyboard,  and  must 
be  a  cataloged  permanent  file.  The  latter  scheme  provides  easy  correction, 
modification,  and  editing  of  the  input  data. 

Each  input  query  is  as  self-explanatory  as  possible.  T!  first  query  at 

the  very  beginning  of  INTERACT I VE-CRKGRO  is  whether  the  user  ined  stress 

intensity  factor  solution  option  is  elected.  If  such  an  opts  is  to  be 
used,  user-defined  stress  intensity  factor  solution  input  is  ired,  first 
for  the  part- through -crack  (if  any)  and  then  for  the  through  lickness 

crack  solution.  A  play-by-play  description  of  this  option 

o  WELCONE  TO  INTERACTIVE-CRKGRO  PROGRAM  VERSION  1.0. 

o  ENTER  FILE  NAME  FROM  CRACK  GROWN  INPUT. 

(Enter  "INPUT."  if  crack-growth  input  data  are  to  be  entered 
from  terminal  keyboard. 

Enter  the  logical  file  name  of  the  external  data  file  containing 
the  crack  growth  input  data. 

NOTE;  All  file  names  must  end  with  a  period;  otherwise,  there  will 
be  a  control  card  error.) 

o  DO  YOU  WISH  TO  ENTER  YOUR  OiVN  STRESS  INTENSITY  FACTOR  EQUATION? 

ENTER  "iTS  OR  "NO."  (Self-explanatory) 

o  IS  YOUR  SIRESS  INTENSITY  FACTOR  FOR  THE  PART-THROUGH  CRACK?  ENTER 
"YES"  OR  "NO".  (Self-explanatory) 

o  ENTER  YOUR  PART-THROUGH  CRACK  STRESS  INTENSITY  FACTOR  EQUATIOTJ 
1 --- . 7- . . . . . 


S 


(If  the  previous  entry  is  YES,  then  this  query  appears.  Enter  the 
user-defined  stress  intensity  factor  equations  in  FORTRAN  statements 
in  Column  7  through  72.  Column  allocations  are  shown  in  numeral  dotted 
points) . 

o  IS  THE  STRESS  INTENSITY  FACTOR  FOR  THE  THROUGH  CRACK?  ENTER  "YES”  OR 
"NO?" 

(Self-explanatory) 

o  ENTER  THE  THROUGH-THE-THICKNESS  CRACK  STRESS  INTENSITY  FACTOR 
EQUATION 

1 - 7 - - - - - 

(Same  as  previous  query) 

To  input  a  user- defined  stress  intensity  factor  solution,  the  user  is 
restricted  to  the  following  variables; 

AOC  Crack  depth  over  crack  length  ratio 

.‘\DT  Crack  depth  over  thickness  ratio 

.\PI.  Crack  depth 

•APREP  Stress  intensity  factor  solution  in  depth  direction/o- 
CPL  Crack  length 

CPREP  Stress  intensity  factor  solution  in  length  direction/^ 

HAI.FW  Half  width  of  structure 
PI  ir  (3.14159) 

Q  Shape  factor 

RADIUS  Radius  of  hole 
T  Thickness 

Since  the  compilation  of  a  new  stress  intensity  factor  subroutine  uses  a 
CX  compiler,  the  user  is  allowed  to  use  any  mathematical  routines  available 
on  the  compiler. 
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INTERACTIVE-CRKGRO  provides  both  one-dijnension  (1-D)  and  two-dimension 
(2-D)  crack  growth  options  for  a  part-through  crack.  For  the  2-D  option, 
when  a  user-defined  K-factor  solution  mode  is  executed,  both  APREP  and 
Q’REP  should  be  defined.  For  the  1-D  option,  only  APREP  needed  to  be 
defined.  In  addition,  a  through-the-thcikness  crack  user-defined  stress 
intensity  solution  is  recommended.  In  the  case  where  a  breakthrough 
(transition  from  part-through-crack  to  through- the-thickness  crack)  has 
occurred,  the  program  will  automatically  switch  from  the  part -through-crack 
subroutine  to  the  through-the-thickness  crack  subroutine.  If  CPREP  is  not 
defined  in  the  through-the-thickness  crack  routine,  QPREP  will  be  set  equal 
to  zero  as  a  result,  there  will  be  zero  crack  growth  in  the  length  direction 
once  the  cracks  have  broken  through. 

Crack-growth  life  conputation  input  consists  of  eight  segments  or 
options  with  labeled  keywords  as  follows: 

1.  TITLE:  Problem  title  and  description 

2.  MATERIAL:  Material  properties  and  crack-growth  rate  constants  and 

parameters 

3.  GEOMETRY:  Structure  geometry  and  load  interaction  effects 

4.  SPECTRLM:  Spectrum  inputs  and  spectrum  information 

5.  CUTPUT:  Establishes  level  of  print  and  graph  desired 

6.  REPLOT:  To  obtain  desired  plots  from  crack-growth  data  saved 

from  previous  analysis 

7.  END  DATA:  Defines  end  of  input  data  and  begins  crack  growth  life 

compulation 

8.  ENDJOB:  Defines  end  of  job  and  terminates  program 

With  the  exception  of  segments  7  and  8,  the  order  of  these  input 
segments  is  inmaterial  to  the  program.  These  segment  inputs  can  be 
overwritten  by  reentering  to  the  segment  where  changes  are  to  be  made  and 
entering  the  necessary  in  puts.  Program  execution  will  be  based  upon  the  last 
input  of  each  segment. 

The  following  are  the  play-by-play  descriptions  of  each  input 
segment  for  performing  crack-growth  life  analysis: 
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o  ENTER  TITLE  CARD  (MAX  OF  70  CHARACTERS) . 

?  Any  alphanumeric  information  which  the  user  desires  to  input  for 
problem  identification. 

2.  OPTION:  MATERIAL 

o  ENTER  MATERIAL  TITLE  CARD  (MAX  OF  60  CHARACTERS) . 

(Self-explanatory) 

o  ENTER  DC/DN  CRACK-GROWTH  RATE  EQUATION  COEFFICIENT. 

?  (Enter  the  crack-growth  rate  equation  coit5)onent  C^  used  in  dc/dN 
equation. ) 

o  ENTER  DC/DN  CRACK- GROOTH  RATE  BQUATICW  EXPONENT. 

?  (Enter  crack-growth  rate  equation  N^  used  in  dc/dN  equation.) 

0  ENTER  DA/DN  CRACK-GROWTH  RATE  EQUATION  COEFFICIENT. 

?  (Filter  crack-growth  rate  equation  coefficient  C^  used  da/dN 
equation) 

o  ENTER  DA/DN  CRACK- GROWTH  RATE  BQUATI(»J  EXPONENTS. 

?  (Enter  crack-growth  rate  equation  exponent  Ng  used  in  da/dN 
equation. ) 

o  ENTTER  CRACK- GROWTH  RATE  EQUATION  EXPWIENTS: 

'M'  for  positive  stress  ratio 

?  (Enter  m,  for  positive  stress  ration  collapsing  factor  m  used  in 
nxxlified  Walker  equation.) 

'Q'  for  negative  stress  ratio 

?  (Enter  negative  stress  ratio  acceleration  index  q  used  in  Chang 
equation.) 
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o  ENTER  THE  FRACTURE  TOUGHNESS  FOR  HALF  CRACK  LENGTH  'C 

?  Plane  stress  or  transition  mode  fracture  toughness,  K  (KSi\/Tn) 
used  for  through-crack  instability  criteria  on  crack  length  C. 

o  ENTER  THE  FRACTURE  TOUGHNESS  FOR  CRACK  DEPTH  ’A' . 

?  Plane  strain  fracture  toughness,  KJ^(ksi^/in.)  used  for  part-through- 
crack  instability  criteria  on  crack  depth  a. 

o  ENTER  mTERIAL  YIELD  STRENGTH  IN  KSI, 

(Self-explanatory) 

o  EI^TER  VALUE  OF  DELTA  K  THRESmiD. 

?  The  threshold  of  Ak  (ksi's/uT. ) 
o  ENTER  THRESHOLD  CONSTANT  A. 

?  (Enter  "ZFJ^”  for  a  constant  threshold,  and  a  non-zero  threshold 
constant  if  a  variable  threshold  is  used.) 

0  ENTER  POSITIVE  STRESS  RATIO  CUTOFF  VALUE. 

?  (Enter  cutoff  value  of  positive  stress  ratio,  Rqu^,  below  which 
material  is  assumed  to  have  no  stress  ratio  layering  effect  on 
crack  growth;  value  between  .0  and  1.) 

o  ENTER  NEGATIVE  STRESS  RATIO  CUTOFF  VAUJE. 

?  (Enter  cutoff  value  of  negative  stress  ratio,  above  which 

material  is  assumed  to  have  no  acceleration  effect  on  crack  growth; 
value  between  -1.  and  0.) 

OPTION:  GE0^ETRY 

o  LOAD  INTERACTION  CONTROL  -  ENTER  ’’NO",  ’TES",  or  "BOIH". 

?  (If  "NO"  is  input,  no  load  interaction  will  be  considered,  and 
compressive  stress  will  be  set  to  zero  in  the  analysis.  If  "YES" 
is  input,  retardation  and  compression  acceleration  will  be 
considered.  If  "BOTH"  is  input, both  analyses  will  be  performed.) 
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o  ENTER  OVERLOAD  SHUTOFF  RATIO. 

?  (Self-explanatory) 
o  ENTER  CRACK  CODE. 

?  (See  Figure  1.) 
o  ENTER  WIDTH  OF  SPECIMEN  (INCH) 

?  (Self-explanatory) 
o  ENTER  THICKNESS  OF  SPECIMEN  (INCH) 

?  (Self-explanatory) 

0  ENTER  RADIUS  OF  HOLE  (INCH,  IF  NO  HOLE,  ENTER  0.0). 

?  (Self-explanatory) 

o  ENTER  TRANSITION  CODE  -  (0)  WHEN  CRACK  EQUALS  THICKNESS 

(1)  BEFORE  CRACK  EC^ALS  THICKNESS 

?  (Control  for  transition  from  part -through  crack  to  a  through  crack) 
0  ENTER  INITIAL  CRACK  LENGTH  (INCH) 

?  (Self-explanatory) 

o  ENTER  FINAl,  CRACK  LENGTH  (INCH),  IF  NO  FINAL  CRACK  LENGTH,  ENTER  0.0.) 
?  (Self-explanatory) 

o  ENTER  NUMBER  OF  LIMIT  STRESSES  (MAXINIJM  OF  7) . 

?  (Self-explanatory) 
o  ENTER  LIMIT  STRESSES  (ksi). 

?  (Self-explanatory) 

o  CONTROL  FOR  DETERMINING  WHEN  INSTABILITY  IS  REACHED  IS: 

1  -  BY  LIMIT  STRESS  INTENSITY  FACTOR 

2  -  BY  MAXIMJM  STRESS  INTENSITY  FACTOR 

(Self-explanatory) 
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OPTION:  SPECTRUM 


o  ENTER  SPECTRUM  TITLE  (MAXIMW  OF  70  CHARACTERS) 

?  (Self-explanatory) 
o  ENTER  SPECTRUM  FACTOR. 

?  (Factor  used  to  scale  the  spectrum,  i.e,,  all  stresses  are 
increased  by  this  factor) 

o  ENTER  RANGE  -  PAIR  COUNTING  CONTTROL. 

0  -  NO  RANGE  PAIR  COUNTING 
1  -  RANGE  PAIR  COUNTING 

?  (Self-explanatory) 

o  ENTER  UNIT  NUMBER  FOR  INPUT  OF  SPECTRUM. 

?  (Enter  ”5”  if  spectrum  is  to  be  input  from  terminal  keyboard. 
Enter  file  unit  number  vhere  spectrum  is  to  be  read  in  from 
external  file.) 

0  THE  FLICHT  SPECTRUM  SEGMENT  TYPES  ARE: 

MAX -MIN  R- DELTA  MFAN  END 

ENTER  SEGMENT  TYPE 

("MAX-MIN"  -  Maximum  and  minimum  stresses  spectrum  are  input 
"R-DELTA"  -  Stress  ratio  and  delta  stress  spectrum  are  input 
"MEAN"  -  Mean  and  alternate  stresses 

"END"  -  Defines  end  flight  segment  spectrum) 

o  ENTER  MISSION-MIX  DEFINITIOJ. 

?  MIS.ST0N  =  NBLKS  (factor  i  and  segment  i) 

NBLKS  -  Nunber  of  times  the  complete  mission  string  is  repeated 
FACTOR!  ■  Number  of  times  individual  flight  segment  is  repeated 
SECMENTi  -  Mission  segment  number  preceded  by  the  letter  "s"  for 
segment 

Sanple:  Mission  *  100  (3S1  +  4S2) 

o  ENTER  NUMBER  OF  FLIGHTS  PER  BLOCK. 

?  (.Self-explanatory) 
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OFriON:  OUTPUT 


o  IS  THE  SPECTRUM  TO  BE  PRINTED?  '’YES"  or  "NO." 

?  (Self-explanatory) 

o  PRINT  CRACK  GROWTH  HISTORY  DATA  IN  INCRE^©f^S  OF - BLOCKS, 

?  Control  for  printing  the  crack  growth  history  in  increments  of 
NB  number  of  blocks. 

0  IS  PLOTTING  TO  BE  DONE?  "YES"  or  "NO". 

?  (If  "NO"  is  input,  the  program  will  skip  the  following  plotting 
questions) 

o  SPECIFY  TYPE  OF  PLOTS  YOU  WOULD  LIKE  -  A  VS  F. 

?  (Enter  "YES"  or  "NO") 

o  DA/DF  vs  F. 

?  (Enter  "YES"  or  "NO") 

0  DA/DF  vs  A. 

?  (Enter  "YES"  or  "NO") 

o  DA/DF  vs  K  /Fit. 

max 

?  (Enter  "YES"  or  "NO") 
o  SPECIFY  TYPE  OF  GRID  DESIRED, 
o  FOR  X-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  F-axis  of  A  vs  F  plot 
o  FOR  Y-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  A-axis  of  A  vs  F  plot 
o  FOR  X-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  F-axis  of  DA/DF  vs  F  plot 
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0  FOR  Y-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  DA/DF  axis  of  DA/DF  vs  F  plot 

o  FOR  X-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  A-axis  of  DA/DF  vs  A  plot 

0  FOR  X-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

?  K  /Fit  axis  of  DA/DF  va  K  /Fit  plot 
max  max 

o  FOR  Y-AXIS,  ENTER  0  -  LINEAR  OR  1  -  LOG. 

DA/DF  axis  of  DA/DF  vs  K „^Flt  plot 
o  SHALL  THE  SLOW  GROWTH  RATE  STEPS  BE  BYPASSED?  "YES"  or  "NO". 

?  (If  "YES"  is  input,  the  program  will  bypass  steps  with  growth  rate 

less  than  10-®  and  AK<2  M  .) 

max  ’ 

0  SHALL  THE  STRESS  INTENSITY  EQUATIONS  BE  PRINTED?  "YES"  or  "NO". 

?  (This  option  does  not  work  if  user-defined  stress  intensity  factor 
solution  is  inputed.  This  option  is  not  available  vrtien  the  user- 
defined  stress  intensity  factor  solution  option  is  invoked.) 

o  SHALL  THE  GROOTH  DATA  OF  THE  FIRST  BLOCK  BE  PRINTED?  "YES"  or  "NO". 

?  (Self-explanatory) 

o  ENTER  NUMBER  OF  STEPS  TO  BE  PRINTED  IN  REQUESTED  PRINT  BLOCK. 

?  (Control  for  printing  intermediate  growth  history  -  i  -  growth  of 
the  first  ith  steps  of  each  "NB"  -  th  blocks  will  be  printed. 

-i  -  growth  of  the  first  i^h  steps  of  each  "NB"  block  will  be 
printed. 

+i  -  growth  of  the  last  i^^^  steps  of  each  "NB"  th  block  will  be 
printed 

0  -  all  steps  of  "NB"  th  block  will  be  printed.) 

6.  OPTION:  END  DATA 

This  will  define  the  end  of  input  data  and  begin  crack-growth  life 
compulation. 


j 
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WHEN  A'?’  IS  PRINTED  ON  SCREEN,  ENTER  A  "P",  AND  PRESS  RETURN,  TO 
START  PRINTING  OR  TO  CONTINUE  PRINTING. 

The  output  consists  of  input  echo  data  and  confuted  results.  The  latter 
can  be  classified  as  either  graphical  displays  or  tabular  displays.  Various 
levels  of  displays  are  provided  for  growth  history  data  \diich  exhibit  a  large 
quantity  of  information.  By  selectively  choosing  the  option,  lisers  can 
create  the  desired  level  of  output  display  vdiich  represents  the  conputed 
results. 

INTERACT I VE-CRKGRO  output  can  be  characterized  into  the  following 
sequential  categories: 

1.  Print  the  input  echo 

2.  Working  spectrum  (optional) 

3.  Conputed  critical  crack  length  in  tabular  form  for  various 
iterations 

4.  Ccmputed  crack-growth  information  in  tabular  form  at  the  desired 
block  intervals 

5.  Crack- growth  summary  table  which  sunmarizes  the  preceding  growth 
history 

6.  Graphical  output  display  of  A  vs  F,  da/dF  vs  F,  da/dF  vs  A,  and 
da/dA  vs  Kmax/F  (optional) 
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Section  IV 


EXAMPLE  PROBLEM 


The  following  exanple  is  selected  to  demonstrate  the  capabilities  of  the 
program.  In  order  to  illustrate  the  procedure  of  executing  the  user-defined 
stress  intensity  factor  solution  option,  the  example  is  to  predict  the  growth 
behavior  of  a  comer  crack  emanating  at  a  fastener  hole  contained  in  a  2219- 
T851  aluminum  plate  subjected  to  transport  spectrum  loading  with  each  flight 
containing  only  the  Ground- Air- Ground  (G-A-G)  cycle  as  shown  in  Figure  2.  The 
stress  intensity  factor  solutions  for  the  comer  crack  at  fastener  hole  dev¬ 
eloped  recently  by  Nevman  (5)  are  employed  as  the  user-defined  K-solutions. 

The  angle  0  was  assumed  to  be  equal  to  90®  in  the  exan^le.  The  following 
is  the  general  stress  intensity  factor  solution  for  comer  crack (s)  at  a 
fastener  hole  as  shown  in  Figure  3,  proposed  by  Newman: 


■'U)"  hw* 


where 

|^Mj+M2  (a/tj^+MjCa/t) 

for 

a/c 

s  1: 

II 

1.13  -  0.09  (a/c) 

o 

M2= 

0.089/(0.2+a/c)  -  0.54 

M3= 

0.5  -  1/(0. 65  +  a/c)  + 

14  (1-a/c)^^ 

g  = 

1  +  jo.l  +  0.35  (a/t)^| 

(1-sin^  )^ 

II 

|(a/c)^  Cos^^  +  Sin^^j 

V4 

Q  = 

1  +  1.464  (a/c)^’^^ 

For  a/c  >  1: 

=  \/t/a  |l  +  0.04  (c/a)  I 

=  0.2  (c/a)^ 

=  -0.il  (c/a)^ 

g  =  1+  [o.l  +  0.35  (q/a)(a/t)^]  (1  -  sinf)^ 
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IP""" 


a 


*^max  =  23.04  ksi 
*^111  i n  =  -11.52  ksi 


Material:  2219-T851  aluminum  plate 


Figure  2.  A  2219-T8S1  Aluminion  Plate  Containing  a  Fastener  Hole 
With  a  Comer  Crack  Subjected  to 
Transport  G-A-G  Cycle  Loadings 
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f0= 


2  2  2 
(c/a)  sin  0  +  Cos  0 


F  =M 
2  4 


Q  =  1  +  1.4.64  (c/a) 

2 


1.65 


1/4 


1  +  O.1(1-COS0) 


0.8  +  0.2  (a/t) 


1/4 


for  a/c  <1:  =  1  +  0.04  (a/c) 

for  a/c  >1:  =  1.13  -  0.09  (c/a) 

1  -  0.15  X^-  3.4ft  ~  4.47  5.52 

^3"  1  +0.13  \2 


^4  = 


X  =  1/  (1  +  c/r  cos  (0.850  )  ) 

IT  D  +  be 


sec 


(  2^  )  sec  (-  2 


W-2c+bc 


y  a/t) 


1/2 


b  =  1 
b  =  2 
1  +  (2-b) 


for  one  comer  crack 

for  two  syimetrical  comer  cracks 


2  + 


ir 

4 


)  (-C. ) 

^  t  ^  R  ^ 


Since  INTERACTIVE  CRKGRO  is  a  two-dimensional  crack-growth  program  for  a 
comer  crack,  the  stress  intensity  factor  solutions  at  the  maximum  depth,  A, 
Ki^,  and  at  the  maximum  length,  B,  K^n,  are  needed  to  be  coded  in.  The  follow¬ 
ing  are  the  inputs  to  INTERACTIVE  CRKGRO: 


UCLCOHE  TO  CRKQRO  ZNTERACTXUe  GRAPHICS  PROQRAH  UERSZON  l.S 


ENTER  FILE  NANE  FRON  WHICH  CRACK  GROWTH  DATA  ARE  INPUT 
FILE  NANE  NUST  DE  FOLLOW  tV  A 
7  INPUT* 

DO  YOU  WISH  TO  ENTER  YOUR  OWN  STRESS  INTENSITY  FACTOR  EQUATION 
ENTER  'YES'  OR  'NO' 

7  YES 

IS  THE  STRESS  INTENSITY  FACTOR  EQUATION  FOR  PTC  OR  TC 
ENTER  'PTC'  FOR  PART-THROUGH  CRACK 
'TC'  FOR  THROUGH  CRACK 

7  PTC 
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ENTER  VOUR  PART-THROUGH  CRACK  STRESS  INTENSITY  FACTOR  EQUATION 


ANGL-Pl/S.  • 

B"1 . 

Ani-i.l3-Q.09«A0C 

Ana — • .  54>G .  89/ ( • .  a4>A0C  ) 

An3-e.5-t./(  .654A0C)«14.t(l.-A0C»ta4. 

G>1.4(0.1'»>.35XA0TSXa.  )XC1.-SIN(  ANGL)  )XXa. 

FANGL- ( AOCXXa . tCOS ( ANGL )«sa. ♦SIN( ANGL  . )XXa . 

0-1.^1.464«A0C««l.fiS 

F 1  -  ( Am  ♦Anas  AOTS  ta .  ♦An3tA0TSS4 .  )SQSFANGL 

An4-i .♦e.e4SAOC 

ANS- 1.  ♦8.  ltd. -COS(  ANGL  nta. 

Ane-Q.a^G.asAOTtte.as 

Fa-An4tAN5tAnfi 

COR-CPL/RAOIUS 

D-a.SRADlUS 

U-a.SHALFU 

ALAnOA- 1 . / ( 1 . ♦CORtCOS( .SSSANGL  )  ) 

An? - 1 . -0 . ISSALAnOA^a . 4eSALAnOAt ta . -4 . 4?tALAnDAtt3 . 43 . sat ALAnDASS4 

F3-An7/d.^e.l3tALAnDAtta.  ) 

An8-i/cos(PitD/(a.tu)) 

An9-l/C0S( (Pl/a )t( (D4BtCPL)/(U-a.tCPL4BtCPL ) >tSQRT( AOT) ) 
F4«SQRT(An8tAn9) 

AniG<a.4PItA0TtC0R/4* 

Anu-a.+pitAOTtcoR/a. 

Ama«soRT(Anie/Anii ) 

Fs<i .♦(a.-B)t(Ania-i . > 

APREP-SQRT(PItAPL/Q)tFltFatF3tF4tF5 

ANGL-8. 

B-1. 

Ani-i.i3-8.09«A0C 

Ana--0,S4+e.89^(o.a+AOC) 

An3-0.5-l./(  .65-«>A0C)t|4.td.-A0C)tta4. 


G- 1 . ♦ ( t . 1 ♦ . aStAOTtta . >t  Cl . -S I N ( ANGL ) ) t ta . 

FANGL • ( AOC t ta . tCOS ( ANGL ) tta . ♦$ IN ( ANGL ) t ta . ) t t a . 
Q-1.^1.464tA0Cttl.65 

FI - ( Ani^ANatAOTtta . ♦An3tA0Ttt4 . > SGtFANGL 

An4-1.40.G4tA0C 

Ans • 1 . ♦G . 1 t (1 . -COS ( ANGL ) >xa . 

Ans -e . 8ts . at AOT t to . as 

Fa-An4tAnstAn6 

COR-CPL/RADIUS 

D-a.tRADlUS 

U-a.tHALFU 

ALANDA • 1 . / ( 1 . TCORtCOS ( . SStANGL  > ) 

AN?- 1 . . 1 St ALAnDA43 . 46tALAnDAtta . -4 . 47t ALANDASta Sat ALARDAt t 4 


P3«  An?/  ( l .  ,  13«ALAnDA>sa .  ) 

An8*i/cos(Px«D/(a.tu)) 

An9-l/C0S(  (PI/a)t(  (D'^B*CPL)/(U*‘2.tCPL'*-BVCPL)  )tSQRT(AOT)  > 
F4-SQRT(An8«An9) 

Anie-2.+PI*A0T*C0R/4. 

Anil-2.<fPISA0T<C0R/2. 

Ani2-SQRT(AniB/Anil) 

F5-l.^(2.-B)«(Ania-l«  ) 
CPREP>SQRT(Pl«APL/Q)<FltFatF3<F4XFS 

END 


Since  there  is  a  possibility  that  the  comer  crack  will  grow  through  the 
thickness  before  the  reaches  its  critical  value,  it  is  necessary  to  input 

the  stress  intensity  factor  solution  for  the  through-crack  counterpart  as 
fol lows ; 


ENTER  VOUR  THROUGH  CRACK  STRESS  INTENSITY  FACTOR  EQUATION 


— -t- --i  — — -  —  J 

B-1.  • 

.  .  U-a.<HALFU 

COR-CPL/RADIUS 

El-l./(l,+COR) 

Ea«t.-0.l5«El«3.4S«El«ta.-4.47«El»t3.«3.5atElS«4. 
F7-Ea/(i+e.i3*Ei»*a. ) 

E3-1./C0S(PI«RADIUS/U) 

E4«B.5tPIt(a.XRADIUS«B«CPL)/(U>a.«CPL'»B«CPL) 

F8«E3<i./C0S(E4) 

E5*a.4B.aS«PltC0R 

E€*a.4e.5«PI«C0R 

e7-SQRT(E5/E6) 

F9-l.  +  (a.-B)»<E7-l,  ) 

CPREP-SQRT(PItCPL)«F7tF8»F9 

END 

lo  execute  the  crack-growth  analysis  portion  of  the  program,  the  follow¬ 
ing  inputs  were  coded  in: 


C  =  3.00(1  X  10* 
a 

C,  =  5.066  X 
b 

10*1° 

n  =  5.8.1 
a 

n,  =  3.  83 
b 

m  =  O.b 
a 

*%  =  0.6 

(j  =0.5 
'a 

qj^  =0.3 

K ,  =  45  ksi  v/  in. 

a 

K  =65  ksi 
"b 

v^in. 

=  48  ksi 

R  =  3.0 
so 

AK^t  =  l.S  ksi  \[Tn. 
th 


A  =  1 


o 

r""  =  0.73 

cut 


R  =  0.99 
cut 
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=  23.04  ksi 


B  =  3.00  in. 


lim 

t  =  0.25  in. 


R  =  0.25  in. 


a.  =  0.05  in. 
1 


a  =23.04  ksi 
max 


C.  =  0.05  in. 

1 

O’  .  =  -11.52  ksi 

min 


In  order  to  provide  a  conplete  picture,  hard  copies  of  the  conputer- 
prompted  statements  and  the  typed-in  inputs  during  the  input  stage  for  exe¬ 
cuting  the  program  to  perform  this  example  analysis  are  presented  in  the 
following: 


INPUT  OPTIONS  AREt 

TITLE  PIATERIAL  GEOMETRV  SPECTRUM 

OUTPUT  REPLOT  END  DATA  ENDJOB 

YOUR  OPTION  IS?  TITLE 
ENTER  TITLE  CARD,  (MAX  OF  70  CHARACTERS) 

?  NEUMANNS  EQUATION  FOR  CORNER  CRACK  FROM  HOLE 
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INPUT  OPTIONS  ARES 

TITLE  MATERIAL  GEOMETRY  SPECTRUM 

OUTPUT  REPLOT  END  DATA  ENDJOB 
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SPECIFY  THE  TYPE  OF  PLOTS  YOU  UOULD  LIKE 
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ITie  outputs  of  the  analysis  are  shown  on  the  following  pages.  The  outputs 
start  with  the  input-echo  and  the  estimation  of  the  critical  crack  length. 
Printouts  of  the  detailed  calculation  in  block  1,  and  so  on,  then  followed. 

A  summary  table  of  the  crack  growth  for  every  block  (21  cycles)  is  shown  as 
the  last  tabulated  output.  For  illustration,  the  CRT  plots  of  crack  length 
versus  number  of  flights  is  also  presented.  Since  one  flight  only  contains 
one  G-A-G  cycle,  the  total  number  of  flights  and  the  total  number  of  cycles 
are  identical. 


C  R  K  G  R  0 

DETAILED  FATIGUE  CRACK  GROUTH  ANALYSIS  PROGRAM 

NEUMAN'S  SIF  EQUATION  FOR  CORNER  CRACK  FROM  HOLE 

CRACK  CODE!  1090  —  USER'S  OUN  PTC  SIF  EQUATION 

LOAD  INTERACTION t  UILLENBORG-CHANG 

DAMAGE  ACCUMULATION:  UROMAN  (LINEAR  APPRX.) 

CRACK  GROUTH  RATE  EG.:  MODIFIED  UALKER 

ANALYSIS  IS  PREFORMED  UITH  A  SPECTRUM  THAT  HAS  BEEN  RANGED  PAIRED 
INSTABILITY  UILL  BE  BASED  ON  MAXIMUM  STRESS 

MATERIAL:  2ai9-T851  ALUMINUM 

DEPTH  LENGTH 

DIRECTION  DIRECTION 

FRACTURE  TOUGHNESS  -45.000  65.000 

GROUTH  RATE  EQ.  CONST.  C  S.0660E-10  5.0660E-10 

GROUTH  RATE  EO.  EXP.  N  3.830  3.830 

GROUTH  RATE  EQ.  EXP.  M  •  .600 

GROUTH  RATE  EX.  EXP.  Q  -  1.000 

YIELD  STRENGTH  •  48.000  KSI. 

DELTH  KTH  -  (1  -  1.0  *  ABS(R))  *  2.500 

+R  CUT-OFF  •  .750 

-R  CUT-OFF  •  -.990 

OUERLOAD  SHUT-OFF  RATIO  •  3.000 

HALF  PLATE  UIDTH  (B)  •  3.000  IN. 

PLATE  THICKNESS  (T)  •  .250  IN. 

RADIUS  OF  HOLE  (R)  ■  .250  IN. 

INITIAL  HALF  CRACK  LENGTH  (C)  •  .050  IN. 

INITIAL  CRACK  DEPTH  (A)  •  .050  IN. 

A/C  RATIO  •  1.000 


ETmLED  FATIGUE  CRACK  GROWTH  ANALVSIS  PROGRAM 
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I I  LED  FATIGUE  CRACK  GROWTH  ANALVSIS  PROGRAM 
NEUMAN'S  SIF  EQUATION  FOR  CORNER  CRACK  FROM  HOLE 
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APPROXIMATE  CRITICAL  CRACK  (A  OR  C )  .  »S38  IN, 

WHEN  K-LIMIT  IS  UITHIN  A  0.5X  TOLERANCE  OF  K  CRITICAL 
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DETAILED  FATIGUE  CRACK  GROWTH  ANALYSIS  PROGRAM 
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